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Abstract 

Wc investigate the properties of cosmological distances in locally inhomogeneous 
universes with pressureless matter and dark energy (quintessence), with constant equa- 
tion of state, px = wxpx, — 1 < wx < 0. We give exact solutions for angular diameter 
distances in the empty beam approximation. In this hypothesis, the distance-redshift 
equation is derived Iron the multiple lens-plane theory. The case of a flat universe 
is considered with particular attention. We show how this general scheme makes dis- 
tances degenerate with respect to wx and the smoothness parameters aj, accounting 
for the homogeneously distributed fraction of energy of the i— components. We analyse 
how this degeneracy influences the critical redshift where the angular diameter distance 
takes its maximum, and put in evidence future prospects for measuring the smoothness 
parameter of the pressureless matter, um- 
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1 Introduction 



Recent measurements of the cosmic microwave background anisotropy ( |de Bernardis et 



al. 2000) strongly support the flat universe predicted by inflationary cosmology. At the 
same time, there are strong evidences, i.e. from dynamical estimates or X-ray and lensing 
observations of clusters of galaxies, that the today pressureless matter parameter of the 
universe ^Im is significantly less than unity {Qm = ^cdm + ^b — [0.30±0.10] + [0.04ib0.01], 
where CDM and B respectively stand for cold dark matter and baryons), while the hot 
matter (photons and neutrinos) is really neglectable. Flatness so requires another positive 
contribute to the energy density of the universe, the so called dark energy - see, for example, 
Silveira & Waga ( | Silveira fc Waga 1997] ) or Turner fc White (| Turner fc White TMTD - 
that, according to the recent data from SNIa ( | Riess et al. 1998 ; Perlmutter et al. 1999 ), 
must have negative pressure to account for an accelerated expanding universe. One of 
the possible solutions of this puzzle is the ACDM universe where the fraction fix of the 
dark energy is supported by a cosmological constant A. A very promising altervative, 
the quintessence, was proposed by Caldwell, Dave & Steinhardt ( paldwell et al. 1998| ) in 
the form of a dynamical, spatially inhomogeneous energy with negative pressure, i.e. the 
energy of a slowly evolving scalar field with positive potential energy. Apart from flatness, a 
scalar field coupled to matter through gravitation can also explain the so called coincidence 
problem, that is, it can provide a mechanism to make the today densities of matter and 
dark energy comparable. As shown in Wang et al. ( Wang at al. 2000 ), a large set of 
indipendent observations agrees with a universe described in terms of pressureless matter 
and quintessence only, a hypothesis that we will make in this paper. 

A scalar field is not an ideal adiabatic fiuid ( | Grishchuk 1994 ; Caldwell et al. 1998) and 
the sound speed in it varies with the wavelength in such a way that high frequency modes 
remain stable still when wx < 0. Moreover, smoothness is gauge dependent, and so a fluc- 
tuating inhomogeneous energy component is naturally defined ( paldwell et al. 1998 ). Inho- 
mogeneities, both in quintessence and CDM, make the relations for the distances derived in 
Friedmann-Lemaitre-Robertson- Walker (FLRW) models not immediately applicable to the 
interpetration of experimental data both in measurements of luminosity distances and an- 
gular diameter distances. The observed universe appears to be homogeneously distributed 
only on large scales (~ 500 Mpc), while the propagation of light is a local phenomena. In 
lack of a really satisfactory exact solution for inhomogeneous universes in the framework 
of General Relativity ( Krasiiiski 1997 ), the usual, very simple framework we shall adopt 
for the study of distances is the on average FLRW universe (Schneider, Ehlers &: Falco 
1992; Seitz, Schneider &; Ehlers 1994), where: i) the relations on a large scale are the same 
of the corresponding FLRW universe; ii) the anisotropic distortion of the bundle of light 
rays contributed by external inhomogeneities (the shear a) is not significant; in) only the 
fraction Oj, the so called smoothness parameter, of the i-component contributes to the Ricci 
focusing TZ, that is to the isotropic focusing of the bundle. The distance recovered in this 
'empty beam approximation', sometimes known as Dyer-Roeder (DR) distance, has been 
long studied (Zel'dovich 1964; Kantowski 1969; Dyer & Roeder 1972, 1973; Linder 1988) and 
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now is becoming established as a very useful tool for the interpretation of experimental data 
(Kantowski 1998; Kantowski & Thomas 2000; Perlmutter et al. 1999; Giovi, Occhionero & 
Amendola 2000). For a very large sample of standard candles, as SNIa for luminosity dis- 
tance, or standard rods, as compact radio sources for angular diameter distances (Gurvits, 
Kellermann & Frey 1999), the effects of overdensity and underdensity along different lines 
of sight between observer and sources balance each other out, and the mean distance of the 
probability distribution of distances for fixed source redshift, corresponding to the FLRW 
distance, can be used (Wang 1999, 2000). Nevertheless, for the present quite poor samples, 
sources are very rarely brightened by the gravitational effect of occurring overdense clumps, 
and so, for propagation of light far from local inhomogeneities, the source appears dimmer 
and smaller with respect to the homogeneous case. So the distance to be used is the mode 
of the distribution, that is, in a very good approximation, the DR distance ( | Schneider et al. 



1992; Kantowski 1998 ). In general, the smoothness parameters for matter and quintessence 



(respectively aju and ax) and the equation of state wx are redshift dependent, but they 
can be considered, in the redshift interval covered by observations, as constant. For the 
dumpiness parameters, the constancy is motivated by the absence of significant variations 
in the development of structures in the observed redshift range. Instead, also in presence 
of redshift variations in wx, only its integral properties have effect on the distance. In fact, 
in flat FLRW models the distance depends on wx only through a triple integral on the 
redshift (Maor, Brustein & Steinhardt 2000). Hereafter, we will consider the Oj and wx as 
constant. For the cosmological constant wx = — 1; a network of light, nonintercommuting 
topological defects ( Vilenkin 1984 ; Spergel &: Pen 1997] ) gives wx = —m/S where m is 



the dimension of the defect: for a string wx = —1/3, for a domain wall wx = —2/3; an 
accelerated universe requires wx < —1/3. 

In this approximation, distances are functions of a family of parameters: and 0,x, 
which describe the energy content of the universe on large scales; wx, which describes 
the equation of state of the quintessence and varies between —1 and 0; two dumpiness 
parameters, aM and ax, representing the degree of homogeneity of the universe and which 
are to be used for phenomena of local propagation. 

This paper is organized as follows: in Sect. 2 we introduce the so called DR equation and 
discuss some features and its solution for quintessence in the form of topological defects; 
Sect. 3 lists the general solution for the DR equation in the flat case and some simple 
expressions in extreme situations; in Sect. 4 the multiple lens-plane theory of gravitational 
lensing is applied for a derivation of the DR equation without using the focusing equation; 
Sect. 5 is devoted to the degeneracy of the distance with respect to the various parameters; 
in Sect. 6 the critical redshift at which the angular diameter distance takes its maximum is 
studied, and, at last, in Sect. 7 we draw the conclusions. 
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2 The beam equation for inhomogeneous quintessence 

In the hypotheses discussed above, the focusing equation (| Sachs 1961| ; | Schneider et al. 



1992) for the angular diameter distance D^^ in terms of an affine parameter A, 

d^DA 



dA2 



\a{X)\'-n{X))DA, 



becomes (see also Linder 1988) 
(PDa 1 



dA2 + t^-^ + [3aMf^M(l + zf + nxax^xil + ^)"^J Da = 0, 



(2) 



where nx = 3{wx + 1) (0 < n^; < 3), and the relation between A and the redshift z, in 
terms of the generalized Hubble parameter 



H{z) = HaJflMil + z)^ + nx{l + zY^ + 0^(1 + z)2, ^^^ = 1-^ 



M 



X, 



IS 



dz 
dX 



• H{z) 



(3) 



(4) 



The isotropic focusing effect in equation is simply represented by the multiplicative 
factor to Da', this coefficient increases with om, olx and nx- 

Equation @ is sometimes called the generalized DR equation. Substituting for A in 
equation (^) by using equation we have 

d^DA 



H\zr- 



dz^ 



■ + 



2H'^(z) IdH'^ 



l + z 



2 dz 



^^ + hl+z) hoM^M + nxax^xil + z^^-^' 
dz 1 ^ 



the initial conditions on equation (|5|) are 

DA{zd,Zd) = 0, 



Da = Q; 
(5) 

(6) 



d_ 

dz 



DA{zd,z] 



1 



z=za '^ + Zd H{zd) ' 

where DA{zd,z) is the angular diameter distance between Zd (that, in general, can be 
different from zero, as occurs in gravitational lensing for the deflector) and the source at z. 
Changing to the expansion factor a = 1/(1 + 2;), equation (||) is 



fDA 
da? 



dDA 
da 



+ 



Da = 0, 



(7) 
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a form which will be useful in the next sections. 

In a generic space-time, the angular diameter distance Da and the luminosity distance 
Dl are related by (Etherington 1933) 

DL = {l + zfDA, (8) 
so that the considerations we will make about Da are easily extended to D^. 
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2.1 Exact solutions for Hk 7^ 

The observational data nowadays available are in agreement with the hypothesis of a flat 
universe, but are also compatible with a non zero, although small, value of ^Ik- A small 
value of is also allowed by the inflationary theory. These circumstances make useful the 
study of the effect of the curvature on the cosmological distances since today technology 
allows to put strong constraints on the cosmological parameters. 
For ax = 1, equation (|^) reduces to 

a^QM + ^Ka + nxa'-^)^ - a (I^m + ^Ra + ^l^xa^-"") ^ + 

cta^ ^ /da 

^aM + ^^^xa3-"^)z)A = 0. (9) 



To solve equation (^), we proceed as in Demiasnki et al. ( Demianski et al. 2000). First, 
we look for a solution in the power form a'^ when J7x = = 0. The parameter s is 
constrained to fulfill the algebraic equation 

- + ^OAf = 0, (10) 

which has the solutions 

s± = ^±^V25 - 24aM = (H) 

When Vix 7^ 0, Q.k 7^ 0, we choose to impose the form Da = cl^ fio) to the solution, being 
/ a generic function. Inserting this expression into equation (P) we have for / 

a [Qm + ^Ka + nxa""-^) ^ " £ (^m{2s - ^) + Qk{2s - l)a (12) 

+ [s{s - 1) - ^s) Qxa'-"^) + {^Ks{s - 2) + ^x{s^ - s{l + ^))a^~^^)f = 0. 

The initial conditions at a = 1 for the auxiliary function / come out from equation (|6|) 
evaluated at 2; = 0, 

/(I) = 0, (13) 



c 

a=l Ho 

Equation (^) is very useful to obtain some exact solutions of the DR equation, corre- 
sponding to integer values of the quintessence parameter nx- Since the solutions for the 
case of the cosmological constant, nx = ( Kantowski 1998 ), and for the pressureless 
matter, nx = 3 ( | Seitz &: Schneider 199^ ), are already known, we will consider only string 
networks, nx = 2 {wx = —1/3) and the domain walls with nx = 1 {wx = —2/3). Let us 
start with the case nx = 2, when equation ( [T^ ) reduces to 

(f f df 
a(ci + C2a)-^ + (c3 + C2C4a)— + C5/ == 0, (14) 



5 



being 

ci = r^M, 

C2 = Oi^ + Ox, 

( 3 

C3 = C2 2s - - 



2 

C4 = 2s — 1, 



C5 



17i^s(s-2) + 17x fsfs-^ + ^)y (15) 



Equation ( |l^ ) is of hypergeometric type, i.e. it has three regular singularities ( |Ince 1956 ), 
and so, for nx = 2, / is the hypergeometric function. If we indicate with fgj^ and fs_ two 
indipendent solutions for, respectively, s = s+ and s = s_, we can write the general solution 
of equation (|9|) for nx = 2 as 

Da = ^+a"*/.+ (o) + A-a'-/.. (o) 

where and A_ are constants determined by the initial conditions. In equation (^) we 
have expressed the scale factor, a, in terms of the redshift. 

Let us consider now the case nx = 1, when the equation for / becomes 

+ (^s {s-2) + (s (s-2) + l)a)f = 0. (19) 



Equat ion ([T7|) is a fuchsian equation with three finite regular points plus a regular singularity 
at 00 ( I Ince 1956 ). The regular points in the finite part of the complex plane are 



01 



02 = , (20) 



a-x = . 

^ 2 



The trasformation y = aja^ sends 02 — > 1 and 0^2 — C = f|- In terms of y, equation ([T7|) is 
V{y-^)[y- C) (^^(2. - -) + ^(2. - l)a22/j 
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which can be reduced to the standard form 



dy'^ \y y-l y-CJ dy 



OXy-q 



where 



7 + -^ 
g (2,-1) a. 



2s 



y{y- i)(y-C) 



/ = o, 



(22) 



7(l + C) + 5C + e, 



tC = — 



(23) 



Equation ( [2^ ) is the Heun equation ( Erdelyi et al. 1955), which is shghtly more comphcated 
then the hypergeometric equation, possessing four points of regular singularity in the entire 
complex plain, rather than three. The constant q is the so called accessory parameter, whose 
presence is due to the fact that a fuchsian equation is not completely determined by the 
position of the singularities and the indices. The Heun equation can be charaterized by a 
V symbol, and the solutions can be expanded in series of hypergeometric functions. Thus, 
the solution of the equation (P) for nx = 1 can be formally written as equation (16), once 
the functions fs^ and fs_ are interpreted as Heun functions. 
In Fig. ??, we plot the solutions found above. 



3 Exact solutions for Q 



K 







The DR equation for a flat universe has already been solved considering the limiting case 
of the cosmological constant ( Kantowski 1998; Kantowski &: Thomas 200C| ; Demianski 
et al. 2000 ). Here, in presence of generic quintessence, we propose the general solution in 
terms of hypergeometric functions and, then, list particular solutions in terms of elementary 
functions. 



3.1 General solution 

When Qk = 0, equation (0) reduces to 



3-nx 



da"^ 



(a^i^M + (1 - ^M)a 
+ i^au^M + ^ax(l - f^Af)a3""^) Da = 0; 



al ^0M + ^(l-f^M)a3~"^ 



dPA 
da 



(24) 
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dividing equation ( p4[ ) by $7^/ and defining /i = ^ ^j'' , we have 

To solve equation (^), we proceed as in Sect. 2.1. First, we look for a solution in the power 
form when ^ = 0. The parameter s is constrained to fulfill equation (|lO|). When 7^ 0, 
we choose to impose the form Da = a^f{a) to the solution, where / is generic. Inserting 
this expression into equation (^5|) and changing to x = a^~"'^ , we have for / 



x(3 - nx)(l + f^^)^+ ((2s - + ^) + (2 - ^nx + 2s^ fix^ ^ 

(s - / = 0. (26) 



nx 



Again, this equation can be solved in terms of hyper geometric functions. Denoting with 
fs^ and fs_ two of such indipendent solutions for, respectively, s = s+ and s = s_, we can 
write the general solution of equation (^4|) as 



Da = A+a''+ fs4x{a)) + A.a'- fU<a)) 



(27) 



(1 + Z)5/4V^' ' \{l + Zf-^x J ' ' V(l + ^)^-"^ 

where and A- are constants determined by the initial conditions. 
3.2 Particular cases 

Once we have the general solution of equation ( p^ in terms of hypergeometric functions, 
we go now to list some expressions of the angular diameter distance in terms of elementary 
functions in two extremal cases. 

3.2.1 Homogeneous universe 

In this case we have that om = = 1, so that the angular diameter distance takes the 
form valid in a FLRW universe, that is 

^ , , c 1 dz , , 

DA{zd,z) = —-— =. (28) 

This is the integral of the differential binomial 

x^(o + 6x^)^ (29) 
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where x = 1 + z, a = OA/, 6 = 1 — $7 a/, ^ = 
this integral in rational form when 

3s - 1 
nx = 



—3/2, h' = nx — 3 and p 
, s eZ-{0}, 



-1/2. We can put 



(30) 



performing the substitutions (Picone & Miranda 1943) 



t = JnM + {i + ^M){i + z)^^ 



when s is even and 



t 



y (l + z)"^-3 

for odd s. Equation ( |30|) includes all and only the rational values of nx for which equation 
( p^ can be solved in terms of elementary functions, nx varies from 2 (wx = —1/3), when 
quintessence evolves like curvature, to 4 {wx = l/3){hot dark matter); for s ±oo, nx 
tends to 3, giving ordinary pressureless matter. For nx = 2 (see also Lima & Alcaniz 2000) 
we get 

DAizd,z) = ^——-^=== Arctanh ( ^^^^^5 ) ; (31) 

Ho (1 + z)Vl - ^AI [ V V 1 - 

we note that with respect to the dynamical equations, a flat universe with nx = 2 behaves 
like an open one with VLk = 1 — r^Af 7^ 0, but, on the other hand, while quintessence 
contributes to the Ricci focusing, a geometric term does not. For nx = 4, it is 



DA{zd,z) 



Ho^m{1 + z) 



a/I + z — ^mz 
(l + z)i/2 



(32) 



Equation ( pj^ holds in the past history of the universe at the epoch of matter-radiation 
equality (zg^ ~ 10^). Other solutions with 2 < nx < 4 are easily found. Even if they can 
be physically interesting when related to other behaviours of the scale factor, they cannot 
explain the today observed accelerated universe. So, we will not mention them here. 

3.2.2 Totally clumpy universes 

We now study very particular models of universe in which both matter and quintessence 
are totally clumped, that is qa/ = ctx = 0. In this case, the DR equation reduces to a first 
order equation and the expression for the angular diameter distance becomes 



c dz 
DA{zd,z) = —{1 + Zd) , 

^0 (1 + z')\ nM{i + z'y^ + nx{i + z T^ 



(33) 



When nx = {wx = —1) and a a/ = 0, the DR equation becomes of the first order 
indipendently of the values of Qx and ax , and so the distance takes the form 



DA{zd, z)=^ii+ Zd) r 

^0 J^dll + z') 



dz 



{l + z)\KlM{l + z)^ + nx 



(34) 
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Once again, in equation ( |3^ ) there is the integral of a differential binomial of the form 
given in equation (|29|), with, this time, a = ^m, b = 1 — ^m, = —7/2, u = nx — 3 and 
p = —1/2. When nx is rational, all and only the solutions of equation (^) in terms of 
elementary functions occur when 



3s — 5 

nx = , s e Z - {0}; 



(35) 



for any such s we can perform the same substitutions already described for homogeneous 
universes in the previous subsection. Now, we have values of nx < 2: for s = 2, 3, 4, 
respectively, we find nx (wx) = 1/2 (-5/6), 4/3 (-5/9), 7/4 (-5/12). For nx = 1/2, it is 



DA{zd,z) 



c 4(1 + Zrf) 
Ho 5(1 - Qm) 



1 - n 



M 



1 - n 



M 



(l+z)5/2 



(36) 



for nx = 4/3, we get 



DA{zd,z) 



Ho 
2 
3 



(1 + Zd) 



M 



1 - J7a/ + 1^m(1 + 2)5/3 (f^^^ _ 1)2 + 



+ -{Qm - if (l - + ^Mil + zf^^) + - (l - + ^Mil + 

and, for nx = 7/4, 



n5/3 



(37) 



Z = Zd 



Da{z) 



Ho 5(1 - n 



M 



1/1-17 



M 



3 \{l + zf/^ 



+ 9. 



M 



3/2 



1 - 9. 



Af1 



M 



(l+z)5/4 



+ 



M 



Other interesting results are obtained when nx = 2 (s = 5) and nx = 4 (s 
nx = 2 (string networks), the angular diameter distance is 

DA{zd,z) = 2^{l + Zd)nli£\{nM-lf + {9m-1?£ 

Ho L 

5 7 



(38) 

(39) 
-5). For 



(40) 



where 



£ = ^l-nM + {l + z)nM; 
for nx = 4 (hot dark matter), it is 



c 1 + Zd 

Ho nljifiM - 1) 



a/1 — f^A^ 

In the limit s — > ±00, nx goes to 3 (CDM). 



+ 



[I + z){i + z - zn 



M 



(41) 
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4 An alternative derivation of the generalized DR equation 



As already shown for a universe with pressureless matter (Schneider & Weiss 1988; Schnei- 
der et al. 1992), it is possible to derive the DR equation from the multiple lens-plane 
theory, without referring to the focusing equation. We want, in the framework of the on av- 
erage FLRW universes, to generalize this result to the case of inhomogeneous quintessence. 
The basic idea is the simulation of the dumpiness by adding to a smooth homogeneous 
background a hypothetical density distribution of zero total mass, which is made of two 
components: a distribution of clumps (both in dust and dark energy) and a uniform neg- 
ative energy density such that the mean density of the sum of both components is zero. 
After such addition, the average properties of the universe on large scales are still that 
corrensponding to the background FLRW model. The gravitational surface density S of 
clumps in a shell of size Az centered on the observer is then 

S = Az^^T^P, (42) 

where the relation between the redshift and the proper distance rprop is that valid in FLRW 
universes, 

dfprop _ c 1 
dz H{z)l + z' ^ ' 

and TPi' is the 0-0 element of the total energy-momentum tensor in clumps, 
rOO(z) = {l-aM)Tl^ + {l-ax)Tf = {I - aM)pM{zy + {I- oix){l + wx)px{zy ■ (44) 

In an on average FLRW universe, the densities of pressureless matter and quintessence are, 
respectively, 

PM{z) = {l + zf^Mpcr, (45) 
px{z) = {l + zT^VLxPcr, 

where pcr = S/Zg / (8vrG) is the today critical density. The dimensionless surface density k 
corresponding to equation (E^) is 



_ 47rG ZJi(z)Di(z,z,) ^ Hi (l + ^)\ ,,^ 

= D,{zs) ''= cHiz) 2 [3(1-«m)1^m+ 



+ nx(l - ax)^x(l + z 



,nx-3 



D,{z)D,{z zs) (46) 
Di{zs) 



where the subscript 1 refers to diameter angular distances in FLRW universes and Zg is a 
hypothetical source redshift. The so constructed spherical shells will act as multiple lens- 
plane. The ray-trace equation which describes successive deflections caused by a series of 
lens planes is ( Schneider et al. 1992| ) 



X; 



Di{zi,z_ 
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where Xi is the bidimensional angular position vector in the ith lens plane and is the 
deflection angle a light ray undergoes if it traverses the ith lens plane at Xj. The lens planes 
are ordered such that Zi < Zj ii i < j. The solid angle distortion is described by the 2x2 
Jacobianes matrices of the mapping equation (|47|) , 



A,.g, (48) 

and by the derivatives of the scaled deflection angle cti = {Di(zi, Zs) / Di{zs))ai, 

U. . p. (49) 
axi 

By taking the derivative of equation ( ^7|) with respect to the indipendent variable xi, which 
represents the angular position of an image on the observer sky, we have the recursion 
relation 

fr[Di{zj)Di{zi,Zs) 

with Ai = X, X being the two-dimensional identity matrix. In our model of a clumpy 
universe, the matrices Ui are given by ( [Schneider fc Weiss 1988 ) 

Ui = -kilAz - %, (51) 

where the first term accounts for the negative convergence caused by the smooth negative 
surface density and % is the matrix that describes the deflection caused by the clumps. In 
the empty beam approximation (light rays propagating far away from clumps and vanishing 
shear), it is 7^ = and then all the Ai are diagonal, Ai = ail. Equation (|50D becomes 



~iDi{zj)Di{zi,Zs) 

where the dependence on Zg drops out in the product of the ratio of distances by ki . In 
the continuum limit, Az —>■ 0, equation (|5^) is 



= 1 + ^'^Xn^'^'i ^iyHy)dy. (53) 
Jo Di{z)Di{y,Zs) 

Multiplying equation ( p2[ ) by Di{z) and letting Da{z) = a{z)Di{z), we obtain, substituting 
for the explicit expression of k given in equation (^) , 



«,.)^ 



3 

-(1 - aM)^M+ 



H{y) 
Di{y,z)DA{z)dy. 
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It is easy to verify that equation is equivalent to the generahzed DR equation equation 
(|5[) with initial condition given for Zd = 0. Changing to Zd the lower limit of the integration 
in equation ( p4[ ) and D^^z) {Di{z)) with Dji{zd,z) {Di{zd, z)), we have the equation for 
generic initial conditions. Equation (^), already derived with a different way of proceeding 
by Linder ( Linder 1988| ), has here been found only using the multiple lens-plane theory. 



Equation ( [5^ ) is a Volterra integral equation of the second kind ( Tricomi 1985| ) whose 
solution is 

DAiz) = Diiz)+ r Hiy,z)Di{y)dy, (54) 
Jo 

where the resolvent kernel H{y, z) is given by the series of iterated kernels 



with 

Ki{y,z) 



Hiy,z)=Y,Ki{y,z), (55) 
1=0 



(?) T7(|f - + ¥(1 - (^x)nxil + 2/)"^^=^] Di{y, z) ii y < z 



elsewhere, 

and the iterated kernels Ki defined by the recurrence formula 

rv 

Ki+i{y,z)= K{y,x)Ki{x,z)dx. (56) 





Since for all i, Ki{y, z) and Di{y) are no negative, we see from equations (|5^)-(|56D that 
the diameter angular distance Dji{z) is a decreasing function of both om and ax, 

Da{z, 4^ ) < Da{z, «1? ) for > «M > (57) 
Da{z,J^'^) < DA{z,aP) for > 



5 Parameter degeneration 

As seen, the consideration of the DR equation in its full generality, with respect to the case 
of a homogeneous cosmological constant, demands the introduction of new parameters. Let 
us study the case of homogeneous dark energy {ax = !)• For ax = 1, equation (^), in 
units of c/Hq, simplyfies to 

Da{z)=Di{z)+ ^i^,(y,z)i?i(y)d?/, (58) 



while equation (|5q) reduces, for y < z, to 



Kiiy,z) = ^{l-aM)nM-^(l + yfDi{y,z). (59) 
2 H{y) 
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Figure 3: The angular diameter distance in ther^M — u^x plane, when aM = oix = 1. The 
distance increases from the top-right to the bottom-left corner, a) We assume z = 0.5; each 
contour is drawn with steps of 0.01. b) We assume z = 1; the step is 0.02. c) We assume 
z = 2; the step is 0.03. d) We assume z = 5; the step is 0.02. The unit of distance is taken 
to be c/Hq. 
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Figure 4: The angular diameter distance in thewx — oiM plane, when = 1 and ax = 1- 
The distance increases from the top-right to the bottom-left corner, a) We assume z = 0.5; 
each contour is drawn with steps of 0.01. h) We assume z = 1; the step is 0.01. It is 
z = 2; the step is 0.02. d) it is z = 5; the step is 0.03. The unit of distance is taken to be 
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Some monotonical properties with respect to the cosmological parameters are then easily 
derived. For an accelerated universe {wx < —1/3), it is 

(9 1 c) 

> 0, j^Diiz) >Oilwx< -1/3, (60) 



dQxH{z) ' d^x 

and so, for every value of the dumpiness parameter ajv/j the angular diameter distance 
increases with increasing Ox, 

d 

Da{z)>Q\Iwx <-l/2,. (61) 



dVLx 



When Wx > —1/3, the inequalities in equation (|60[ ) are reversed and the distance decreases 
with increasing 0,x- With respect to the equation of state wx, it is 

< 0, ^ Di{z) < 0; (62) 



dwx H{z) dwx 

and so 



^ -Da{z) < 0; (63) 



dwx 

large values of the distances correspond to large negative values of the pressure of quintessence 
so that, for fixed VI m, and aM, the angular diameter distance takes its maximum when 
the dark energy is in the form of a cosmological constant. 

We now want to stress the dependences of the angular diameter distance on wx, 
and au in flat universes with ax = 1- As Fig.|| and Fig.|2| show, the angular diameter 
distance is degenerate with respect to different pairs of parameters, since the distance in 
the ACDM model with = 0.3 is not distinguishable, within the current experimental 
accuracy (Perlmutter et al. 1999), from the one in a FLRW universe with less pressureless 
matter but a greater value of wx or from an inhomogeneous universe with greater wx and 
the same content of matter. 

In Fig.|3|, we plot the degenerate values of the distance in the Q-m — wx plane when 
universe is homogeneous for four different source redshifts: as expected, the dependence of 
the distance on the cosmological parameters increases with the redshift of the source. A 
general feature is that the distance is less sensitive to the components of the universe when 
Qm is near unity and wx goes to 0. This is easily explained: when Qm is large, quintessence 
density Ox is not, and the pressureless matter characterizes almost completely the universe; 
moreover, a value of wx near zero describes a dark energy with an equation of state very 
similar to that of the ordinary matter. So, increasing wx mimics a growth in ^m- On the 
other side, for low values of Om (wx) the distance is very sensitive to wx {^m) and this 
effect increases with the redshift. We see from Fig.^ that the effects of wx and Qm are of 
the same order for a large range of redshifts. 

In Fig.^ we compare, for Qm fixed to 0.3 and for different source redshifts, the compelling 
effects of aM and wx on the distance. When um goes away from the usually assumed value 
(ajvf = 1), once fixed the redshift, the distance increases; on the contrary, for wx that 



17 




Figure 5: For flat homogeneous universes Zmis determined by the intercept between the 
angular diameter distance and the always decreasing Hubble distance. The values on the 
ordinate axis are in units of c/Hq. It is l^Af = 0.3, wx = — 1- 



goes away from the value corrensponding to the cosmological constant (wx = — !)> the 
distance decreases. The dependence of the distance on um increases very rapidly with z, 
and, when z = 5, the effects of um and wx are of the same order. From Fig]^ we deduce 
that the dependence on um increases when wx goes to —1, since values of wx near zero 
have the effect to smooth the universe. In fact, when wx = — 1, both a fraction om of 
the pressureless matter and of the cosmological constant are uniformly distributed; when 
wx 0, quintessence behaves like ordinary matter, and so, for the same value of om, 
the pressureless matter homogeneously distributed is om^m + = 1 — (1 — aM)^M- 
Intermediate values of wx interpolate between these two extreme cases. 



6 The critical redshift 

The critical redshift at which the angular diameter distance of an extragalactic source takes 
its maximum value has already been studied for the case of a flat ACDM universe by Krauss 
& Schra mm (|Krauss &: SchrammlQQ^ ) and for a flat universe with quintessence by Lima &: 
Alcaniz ( Lima fc Alcaniz 20001) . In this section, we will find again their results with a new 



approach and will extend the analysis to inhomogeneous flat universes. 

As can be seen cancelling out the derivative of the right hand of equation ( p8[ ) with 
respect to z, the critical redshift Zm for a flat homogeneous universe occurs when 

DA{Zd,Zm) = JT7 (64) 

SO that the angular diameter distance between an observer at z = z^ and a source at z^ 
is equal to the Hubble distance for z = Zm, ats you can see in Fig.^. Equation (|6^ ) is an 
implicit relation that gives the dependence of Zm on z^, and wx ■ Throughout this 
section, we will put z^ = 0. In Fig.^ we show Zm for a homogeneous flat universe. For a 



18 



0.2 0.4 0.6 O.S 1 



Figure 6: Contours of equal on the ^Im — ij^x plane for flat homogeneous universes. Each 
contour is drawn with a step of 0.05. 
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Figure 8: Contours of equal Zm on the Q.m — aMplane for universes with a cosmological 
constant. Contours are drawn with steps of 0.03. 




Figure 9: Contours of equal Zm on the wx — OiM plane for flat universes with J7m = 0.3 and 
ax = 1- Contours are drawn with steps of 0.03. 
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given value of wx (^^m), Zm decreases with increasing Qm {wx)', when Qm = 0, Zm diverges 
for Wx = — 1, but also a small value of JIm is sufficient to ehminate this divergence (see 
Figj^. The minimum value of corresponds to the Einstein-de Sitter universe {^m = 1 
or Wx = 0), when = 1-25. As you can see from Fig.^, for values of wx in the range 
(—1,-0.8), once fixed fijvf, z^ is nearly constant and this trend increases with ^Imj on 
the contrary, for small Om (~ 0.4) and wx ~ —0.4, z^ is very sensitive to wx- The small 
changes of z^ in the region of large 0,m and wx are explained with considerations analogous 
to those already made in the previous section for the values of the distance in the — wx 
plane. 

Let us go now to analyse the effect of aM on Zm- By differentiating equation (^3|) 
and equation (^) with respect to z, we see that the derivatives are zero only for z ^ oo: 
i.e., in flat universes with totally inhomogeneous quintessence or in a generic model with 
cosmological constant, the critical redshift is not finite when a^./ = 0. So with respect to 
Zm, a totally clumpy universe, indipendently of Um and wx, behaves like a FLRW model 
completely dominated by the vacuum energy. In fact, the cosmological constant, differently 
from dark energy with wx > — does not give contribution to the Ricci focusing and the 
same occurs for the pressureless matter with aM = 0. In Fig.^ we show Zm in the — om 
plane for wx fixed to —1. The critical redshift decreases with increasing fijv/ and aM, and 
takes its minimum for the Einstein-de Sitter universe {^m = ctM = 1), that is when the 
focusing is maximum. On the other side, Zm is very sensitive to aM, especially for large 
values of Qm since aM appears in the DR equation as a multiplicative factor of ^m- For 
Qm = 0.3, Zm = 1-61 and 3.23 for, respectively, aM = 1 and 0.2, a variation of 100%. 
So, combining different cosmological tests to constrain the other cosmological parameters, 
we can use the redshift-distance relation to guess the smoothness parameter aM in a quite 
efficient way. 

We conclude this section comparing the influence of aM and wx on the critical redshift. 
Fig]^ displays Zm in the aM — wx plane, for Qm fixed to 0.3 and with ax = 1- As expected, 
Zm increases when the focusing decreases, that is for small values of aM and wx- We can 
see that the effects of aM and are of the same order. 

7 Conclusions 

The question of light propagation in an inhomogeneous universe is an open topic of modern 
cosmology and the necessity of a positive solution grows up with the increasing means 
of experimental technology that begins to explore very high redshifts. In particular, the 
cosmological distances, which provide important probes of the universe, are very sensitive 
to cosmological inhomogeneities. In fact, in the present quite poor data samples, sources 
usually appear dimmer and smaller with respect to homogeneous models for the systematic 
effect of under-densities along the different lines of sight. Here, we have investigated the 
properties of the angular diameter distance in presence of either clumpy pressureless matter 
and inhomogeneous dark energy using the 'empty beam approximation'. 

The equation for the angular diameter distance with respect to the redshift has been 
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found in a way that is indipendent from the focusing equation. The multiple lens-plane 
theory allows to derive the DR equation for inhomogeneous quintessence, in a way that 
makes clear the importance of the 'empty beam approximation' in the gravitational lensing. 

We have given useful forms for the distance. For non flat universe, we have studied 
the case of cosmic string network, when the angular diameter distance is expressed in 
terms of hypergeometric functions, and an accelerated universe with wx = — 2/3 (domain 
walls), when the distance-redshift relation is given in terms of Heun functions. For the very 
interesting case of flat universes with inhomogeneous quintessence, we have obtained the 
solution of the DR equation in terms of hypergeometric functions. Then, we have listed 
the expressions for the distance when it takes the form of elementary functions for two 
particular cases of flat models: strictly homogeneous (we have considered cosmic string 
networks again) and totally inhomogeneous universes, i.e. both pressureless matter and 
dark energy solely in clumps. 

As it could be reasonably expected, the angular diameter distance is degenerate, also 
considering a smooth quintessence, with respect to cosmological parameters: although there 
is a strong dependence both on the smoothness parameter om and the equation of state 
Wx, the intrinsic systematic errors prevent us to univocally determine the cosmological 
parameters only from measurements of distance. Nevertheless, combinig these data with 
indipendent observations like that concerning the cosmic microwave backgroung radiation, 
we can estimate the smoothness parameter ckm by determining the critical redshift, for 
example. 

The effect of wx and aM, especially at high redshifts, are of the same order but very 
different. An increase in wx mimics a larger il-M, while a more inhomogeneous universe 
(small om) mimics the effect of a cosmological constant. This last statement agrees with 
what found by Celerier ( pelerier 2000| ), in a very different context, that is by relaxing the 
hypothesis of large scale homogeneity for the universe. We have obtained a similar result 
by adopting the cosmological principle. 
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